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Abstract 


We treat general relativity as an effective field theory, obtaining 
the full nonanalytic component of the scattering matrix potential to 
one-loop order. The lowest order vertex rules for the resulting effective 
field theory are presented and the one-loop diagrams which yield the 
leading nonrelativistic post-Newtonian and quantum corrections to the 
gravitational scattering amplitude to second order in G are calculated 
in detail. The Fourier transformed amplitudes yield a nonrelativistic 
potential and our result is discussed in relation to previous calculations. 
The definition of a potential is discussed as well and we show how the 
ambiguity of the potential under coordinate changes is resolved. 


1 Introduction 


The idea that a field theory need not be strictly renormalizable in the tra¬ 
ditional sense yet can still yield useful quantum predictions when treated as 
an effective field theory |jl| has been clearly demonstrated in chiral pertur¬ 
bation theory and in other applications]^. Quantum loop calculations lead 
to well defined results in the low energy limit. Interestingly, such methods 
can also be applied to general relativity. As an effective field theory, the 
quantization of general relativity can be carried out in a consistent way, 
since the troublesome singularities which occur for various types of matter 
sources in traditional renormalization schemes ililH can be absorbed 
into phenomenological constants which characterize the effective action of 
the theory. This effective field theory approach offers then a possible way 
around the familiar renormalization difficulties of general relativity in the 
low energy regime and, using this approach with background field quanti¬ 
zation 1^ of general relativity, one of us (^, ]T| some years ago derived the 
leading quantum and classical corrections to the Newtonian potential of two 
large non-relativistic masses. This calculation has since been the focus of 
a number of publications [0, 0, H, |T^, and this work continues, most 
recently in the paper [^. Unfortunately, due to difficulty of the calculation 
and its myriad of tensor indices there has been little agreement among these 
various authors. The classical component of the correction has previously 
been discussed in the papers by|l, 0, [I^, [I^ , p^] , and here there is gen¬ 
eral agreement although, as we shall discuss, there exists an unavoidable 
ambiguity in defining the potential. The basic disagreements lie rather in 
the quantum corrections and in the present paper we shall present what 
we believe to be the definitive result for the leading classical and quantum 
corrections of order G^, using the full scattering amplitude as the definition 
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of the non-relativistic potential. 

We note that, as a prelude to this effort, in a recent paper I'M] two of us 
have dealt with the the quantum and classical corrections to the Reissner- 
Nordstrom and Kerr-Newman metrics of charged scalars and fermions. Such 
quantum and classical corrections have also been considered from the view¬ 
point of a scattering potential in a paper by one of us. Recently we have 
also calculated the full classical and quantum corrections to the Schwarzshild 
and Kerr metrics of scalars and fermions |26| and have shown in detail how 
the higher order gravitational contributions to these metrics emerge from 
loop calculations. In the present paper then we consider the corresponding 
calculation of the full scattering amplitude. 

Of course, treating general relativity as an effective field theory is carried 
out at the cost of introducing a never ending set of additional higher deriva¬ 
tive couplings into the theory. In this sense Einstein’s general relativity is 
still a perfectly valid theory for gravitational interactions—although now it 
represents only the the minimal theory. At some stage additional derivative 
couplings must be appended to the Einstein action, signifying manifesta¬ 
tions of the higher energy component of the effective field theory. However, 
the low energy scattering potential is free from these new couplings and 
represents a model-independent result for quantum gravity. 

This calculation is possible because the post-Newtonian and quantum 
corrections which we consider are determined fully by the non-analytic pieces 
of the one loop amplitude generated by the lowest order Einstein action. Of 
course, in order to deal with the ultraviolet divergences which arise at one 
loop, one must renormalize the parameters of higher derivative terms in 
the action. However, such pieces will only affect the analytic parts of the 
one-loop amplitude, and will not contribute to our potential. 

We will employ the same conventions as in our previous papers, namely 
(h = c = 1) as well as the Minkowski metric convention (-1-1, —1, —1, —1). 
We will begin in section 2 with a very short introduction to the effective 
held theory quantization of general relativity and focus here especially on 
the distinction between non-analytic and analytic contributions to the scat¬ 
tering amplitude. We also include a discussion of the dehnition of the non- 
relativistic potential. 

Next in section 3 we evaluate the diagrams which contribute to the scat¬ 
tering and examine in detail the results for the various components. The 
resulting nonanalytic piece of the scattering amplitude is then used in order 
to construct the leading corrections to the nonrelativistic gravitational po¬ 
tential. We also discuss our result in relation to previous calculations and 
attempt to sort out the various inconsistencies in the published numbers. 
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Finally in a concluding section we summarize our findings. 


2 Review of general relativity as an effective field 
theory 

We begin with a brief review of general relativity, the Lagrangian of which 
(not including a cosmological term) is 



( 1 ) 


where k? = 32ttG is the gravitational coupling, — dpT^a + 

^aa^Zp — is the curvature tensor, and g denotes the determinant 

of the metric field 5 ^,^. Here ^/Imatter is a covariant Lagrangian for the 
matter fields, and in principle any type of matter field could be included. 
This action defines the classical theory of general relativity. 

In order to treat Eq. |I| as an effective field theory one must include 
all possible higher derivative couplings of the fields in the gravitational La¬ 
grangian. In this way any field singularities generated by loop diagrams can 
be associated with some component of the action and hence can be absorbed 
via a simple redefinition of the coupling constants of the theory. Treating 
all such coupling coefficients as experimentally determined quantities, the 
effective field theory is then finite and contains no singularities at any finite 
order of the loop expansion. 

We can consequently write an effective action for pure general relativity 
as 




where the ellipses denote that the effective action is in fact an infinite series— 


at each new loop order additional higher derivative terms must be taken into 


account. This Lagrangian includes all possible higher derivative couplings, 
and every coupling constant in the Lagrangian is considered to be deter¬ 
mined empirically. Similarly one must include higher derivative contribu¬ 
tions to the matter Lagrangian in order to treat this piece of the Lagrangian 
as an effective held theory. Details of such considerations can be found in 
the papers [^, 

In our calculations we will consider only the non-analytic contributions, 
which are generated by the propagation of two or more massless particles 
in the Feynman diagrams. Such nonanalytic effects are long-ranged and. 
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in the low energy limit of the effective field theory, they dominate over the 
analytic contributions which arise from the propagation of massive modes. 
The typical non-analytic terms we will consider are of the type: 1/y—^ 
and In —while the typical analytic contribution is a power series in q. The 
feature that the analytic contributions originate from the propagation of 
massive particles while the non-analytic effect comes from massless prop¬ 
agation can be seen directly by Taylor expanding the two types of prop¬ 
agators. Indeed the massless propagator l/q‘^ cannot be expanded in a 
series, while we have the obvious representation for the massive propagator: 
l/(g^ — m?) = —— q^/m? + ...). Thus the massive parts of the 
diagrams will always be expandable and hence analytic—while the massless 
contributions have the possibility of generating nonanalytic components. It 
should be noted that such nonanalytic pieces of the scattering amplitude are 
essential to the unitarity of the S-matrix. 


2.1 Definition of the potential 


Before proceeding to the actual calculations, it is important to note that the 
precise definition of a potential in a relativistic quantum field theory such 
as general relativity is not obvious. In the original papers, the one-particle- 
irreducible potential was calculated H], |l^. However, additional diagrams 
are required in order to relate this quantity to physical observables. Sub¬ 
sequent work has considered alternative definitions. Clearly a gravitational 
potential should be gauge invariant in order to make physical sense, but one 
can choose between various definitions of the potential depending on the 
physical situation, how the total energy is defined, etc. In QCD one can al¬ 
ways turn to the venerable Wilson loop description but there exists no stan¬ 
dard Wilson loop description for the gravitational potential, although some 
work has been done in this direction using the Arnowitt-Deser-Misner 
formula for the total energy of the gravitational system. Such an approach 


to the gravitational potential has been taken in the paper |11]. For yet an¬ 
other approach to the Newtonian potential, see [^]. An alternative path is 
to use the scattering amplitude itself to define the potential. This descrip¬ 
tion of the potential seems to us the simplest and most intuitive picture and 
has been employed by a number of authors [17, H, 0, Herein 

then we shall also use the full scattering amplitude in order to represent the 
potential, defining 


{f\T\i) = {27,)H^^\p-p'){M{q)) 

= -(27rM(i?-F;0(/|H(q)K) (3) 
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where p,p' is the incoming, outgoing four-momentum. The corresponding 
coordinate space representation can be found by taking the nonrelativistic 
limit and Fourier-transforming, yielding the result 

and will serve as our definition of the nonrelativistic potential. 

It should be noted, however, that this is not the only way in which to 
define the potential F(q) in terms of the scattering amplitude. One could, 
for example, subtract off the second order Born contributions, which would 
lead to the nonrelativistic potential used in bound state quantum mechanics 
and would be equivalent to using the prescription 


i{f\T\i) = 


X 


-27ri5{E - E') 


{f\Vbs{q)\i) +Y1 

n 


{f\Vbs{q.)\n){n\Vbs{ci)\i) 
E - En + ie 



(5) 


The definition of the bound state potential is discussed in detail in |^]. In 
particular, in a Hamiltonian treatment there are also terms in the Hamil¬ 
tonian involving G]P‘jr that contribute at the same order. The relation of 
the bound state potential I4s(q)) in Einstein-Infeld-Hoffmann coordinates, 
to the lowest order scattering potential is 


Vbs{r) = V{r) + -- 3-2 - 


( 6 ) 


2.2 The diagrams contributing to the nonanalytic compo¬ 
nent of the scattering matrix 

We will consider here only the nonanalytic contributions from the one-loop 
diagrams. Since many diagrams yield purely analytic contributions to the S- 
matrix, such diagrams need not be considered and will be omitted from the 
beginning. The diagrams which do yield nonanalytic contributions to the 
S matrix amplitude are those containing two or more massless propagating 
particles. A typical such amplitude will be of the form 

M ^ (a + Bq^ -k ... -k + fdiK^ In(-g^) -k 132^^—^= + . • • ) (7) 

Here the coefficients A,B,... correspond to analytic pieces which are of no 
interest to us, as these terms will only dominate in the high energy regime of 
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the effective theory. Rather, the a, Pi, f32, ■ ■ ■ coefficients correspond to the 
nonlocal, nonanalytic contributions to the amplitude and are the ones which 
we seek. In particular, the Pi,P 2 terms will yield the leading post-Newtonian 
and quantum corrections to the potential. 

3 Results for the Feynman diagrams 

In this section we will present our results. All diagrams have been per¬ 
formed both by hand and by computer. In order to evaluate the diagrams 
by computer, an algorithm for Maple 7 (TM)[| was developed. This program 
contracts the various indices and perform the loop integrations. All results 
obtained this way were confirmed results obtained by hand. The resulting 
amplitudes were then Fourier transformed to produce the scattering poten¬ 
tial, and only the nonanalytic pieces of the amplitude were retained. For 
this part of the calculation, the following Fourier integrals are useful: 



After this brief introduction, we proceed to give the results for each diagram 
in turn. Note that the basic vertices needed for our calculation are given in 
Appendix A. 

3.1 The tree diagram 

The result for this diagram in the nonrelativistic limit is the well-known 
lowest order tree-level result which yields the Newtonian potential. We 
define the diagram using the Feynman rules as: 



where q = ki — k 2 = — k^. Contracting all indices and taking the nonrel¬ 

ativistic limit, we find 



( 10 ) 


^Maple and Maple V are registered trademarks of Waterloo Maple Inc. 
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(mi) (m2) 

^ ki fcsX 


1 (a) 


Figure 1: The tree diagram giving Newtons law. 

whose Fourier transform produces the scattering potential 

Gmim-2 

Viia){r) = -^- 

which is indeed the familiar Newtonian form. 

3.2 The box and crossed box diagrams 


( 11 ) 


< 

k2 ki 

(mi) ' ' (m2) 

/ \ 

(mi)^^(m2) 
k± ^ 

/ \ 

2 (a) 

2 (b) 


Figure 2; The box and crossed box diagrams which contribute to the non- 
analytic component of the potential. 

We can write the contributions of these two diagrams as 

M2a = J + 

X T^^iks, /cs - I, m2)Tl\k-i - I, /c4, m2) (12) 


i 

i 

i'PfiuafS 

i'Ppa-yS 

-{ki + l)"^ — ml- 

-(/c3 — l)'^ — m\- 

[ p \ 
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for the box and 


f dH 

^26 = J ^2^)4 ^ ^2, ?7J,i) 

X r7'^(/c3 ,1 + ki-,'m2)T^^{l + A:4,m 2 ) (13) 


i 


i 


PP iiual3 

i'P pa'yS 

-{ki + /)2 — ml- 


-{1 -I- kY‘^ — ml- 


[ P \ 

Hl + q)^\ 


for the cross box. These diagrams are among the most challenging that 
we will encounter, due to the rather complicated integrals—containing four 
propagators—which must be evaluated. However, those pieces of the am¬ 
plitude which are loop momentum-dependent simplify, due to the feature 
that the external particles are on shell and by the fact that we are only 
seeking the nonanalytic pieces of the scattering amplitude. This allows re¬ 
duction of parts of the amplitude initially having four propagators to pieces 
where effectively only three or two propagators remain. An example of this 
simplification can be seen by the replacement of the integral 


dH __ 

( 27 r )2 P{1 -|- + fei)^ — mf)((/ — /cs)^ — m^) 

_ 1 r JH _ {{I + ki)^) - P - mj _ 

2 J ( 27 r )2 P{1 + qP{{l -h kiP - ml){{l - k-^Y - ml) 


(14) 


by 

1 
2 

because the P part will not contribute to the nonanalytic component. The 
integrals with three or two propagator terms are explicitly given in Appendix 
B. Another simplification arises when the momentum ( 7 ^ contracts with a 
loop momentum P in the numerator. An example is 


dH 


( 27 r )2 p[l q)‘^{{l — ^3)2 — ml) 


(15) 


dH _ Pj _ 

( 27 r )2 P {1 + q)‘^{{l + ki)‘^ — mf)((/ — ks)'^ — ml) 

^ 1 f _ {{I + q)^) - P - qf _ 

2 J ( 27 r )2 P(l + q)‘^{{l -h ki)"^ - mf)((/ - ^3)^ - ml) 


(16) 


which simplifies to 


1 

2 


d'^l q^ 

(27r)2 P{1 + q)‘^{{l — ks)^ — ml) 

8 


(17) 



Via these simplifications one can reduce the box and cross box amplitudes to 
a reduced piece consisting only of integrals with two or three propagators and 
a component with the basic form of the box and crossed box integrals, he., 
with no loop momentum terms in the numerator. Then, using the integrals 
presented in Appendix B, performing the above-described contractions in 
the two diagrams, and taking the nonrelativistic limit we end up with the 
result 

46 

^ 4 "(J+ 2 ( 6 ) log f (18) 

for the momentum-reduced component of the box plus cross box and 


^ 2 ia)% 2 (b) log f (19) 

for the irreducible piece. Fourier transforming, we find then the scattering 
potential contributions 


for the reducible component and 


23 mim2G‘^ 
3 vrr^ 


T/irred 

^2{a)-\-2{b) 


(r) = -8 




( 20 ) 


( 21 ) 


for the irreducible piece so that the total result for the box and cross box 
contribution to the potential is 


T/tOt 

'^2(a)+2(b) 


(r) = 


47 mim2G‘^ 
3 vrr^ 


These results are in agreement with those of 


H- 


( 22 ) 


3.3 The triangle diagrams 

The next pieces we will consider are the set of triangle diagrams, for which 
we find 


M3{a){q) 


dH 

{2-kY 


-I- ki,mi)TY^{I -|- fci,/c2, ? 7 ri)r^^'’''^(/i; 3 , ^4, m2) 


i'POLp'^S 

i'Pfiuap 

i 

l{l + qy\ 

[ P \ 

-{1 + kip — m\- 


( 23 ) 
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Figure 3: The set of triangule diagrams contributing to the scattering po¬ 
tential. 


and 


= j - l,m2)Tf{k3 - l,k4,m2)T^'''^^{ki,k2,mi) 



OL0''y5 

i 

[ P J 


-{1 — ^ 3 )^ — m-i- 


(24) 

The calculation of such diagrams yields no real complications— the integrals 
needed are quite straightforward and are presented in Appendix B. However, 
a significant simplification results from the use of the identity 


{kl,k 2 ,mi) = (^ 1 , ^ 2 , "ll) 

and, taking the nonrelativistic limit, we find for these two pieces 

-2 , 


(25) 


M3(a){^ = -8G^mim2 Qlog^ + 
M 3 {b) Q log f ^ 


(26) 


Our results for these diagrams agree with those of |^, and the Fourier 
transformed result is 

^ G^mim2(mi + m2) 

h3(a)-|-3(fe) (^) 4 2 +28 _ ^ [27) 




3.4 The double-seagull diagram 

We have for the double-seagull term 


M4(a) (q) = 


1 


dH 


2! J (27r)4 


T^d'y^[ki,k2,mi)T2^^''{k3,ki,m2) x 


'^P OLjdpiU 

ip'ySap 

l{l + q)^\ 

[ p \ 


( 28 ) 
























4(a) 


Figure 4: The double-seagull diagram contribution to the scattering poten¬ 
tial. 


The double-seagull loop diagram is quite straightforward, and is simplified 
by use of the identity Eq. 25. Note, however, that there exists a symmetry 


factor of 1/2!. The resulting amplitude is found to be 


= 44G^mim2logg^ 


(29) 


whose Fourier transform yields the double-seagull contribution to the po¬ 
tential 


F4(,)(r) = -22 


mim2G‘^ 

vrr^ 


(30) 


3.5 The vertex correction diagrams 


f^2 V 

(mi) '( {m 2 ) 


5(a) 



5(c) 


k2 f k4 

(mi) ( ("^ 2 ) 

\k3 

/ \ 


5(b) 



Figure 5; The class of the graviton vertex corrections which yield nonana- 
lytic corrections to the potential. 
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There exist two classes of vertex correction diagrams: For the massive 
loop diagrams, shown in Figures 5(a) and 5(b) we have 

= ■^^^'rr^(^i,^2,mi)rf''(/c3,/ca - /, m2)rf'^(/c3 - 


Xk^u 

(pepa 

i'Pa.13'^6 

i 

[ P J 

Ui + q^i 

V g2 J 

-{1 — k^Y — 1172- 


(31) 


and 


M 5 (b)iq) = j {ki,l + A:i,mi)Tf''(Z + A:i,A; 2 ,mi)Ti^'"(/c 3 ,/c 4 ,m 2 )r/^'^^'^^^(-Z,g) 



i'Ppupa 

i'P (peXn 

i 

[ p \ 

Hl + qy\ 

. - 

-{1 -|- fei)^ — mf. 


(32) 

while for the pure graviton loop diagrams shown in Figures 5(c) and 5(d)we 
have 

1 ^ 4 (c)(g) = ^ j ^^j:j4 7 - 2 ^''"'^(fe 3 ,fc 4 ,rn 2 )rf^(A;i,fe 2 ,mi)T 3 ^‘'^'"^^'^\z,-g) 


PP fluXK 

paep 

ip 

[ P \ 

ki + q^i 

. q^ . 


(33) 


and 


dH 




ippu'yS 

ip paal3 

ip Xtiep 

[ p \ 

ki + q?i 

- - 


(34) 

The vertex correction diagrams are certainly the most challenging to per¬ 
form, disregarding the box and crossed box diagrams, and the results go 
back to the original calculation of |^, |I^— however, due to an algebraic er¬ 
ror, the original result quoted for the such diagrams was in error. Since that 
time, the results for the gravitational vertex corrections have been checked at 
length in various publications |T^, |l^; however, until the correct forms 
have not been given. Using the results of [^] and taking the nonrelativistic 
limit, we find the amplitudes, 


M5{a)+5{b)i^ = 2G‘^ 17111712 

52 


7r2(mi-km2) 5 .g 

- + 3 logg 


^5(c)+5(d)(^ = — -G^mim2logq 


\q\ 

2 


(35) 
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whose Fourier transform yields the corrected results for the vertex modifi¬ 
cations to the scattering potential: 


^5{a)+5(b) (^) 


G‘^mim2{mi -|- m 2 ) 


5 mim2G‘^ 
3 vrr^ 


(36) 


and 

26 mim2G^ 

H(c)+5(d) (^’) — 2 ^^3 

where again we note the presence of a symmetry factor 1/2! in the case of 
the pure graviton loop diagram.^ Our results for these diagrams are not in 
agreement with previous calculations. For the vertex correction diagrams 
5(a) and 5(5), there is total disagreement with our result. However, we have 
performed a very detailed analysis of the vertex diagrams in our compan¬ 
ion paper |^, showing how the classical ingredients match in detail those 
required by the classical Schwartzschild metric. The quantum vertex result 
of |l^ is the same as that of iQ and the latter yields an incorrect classical 
metric (Ref. |l^ does not display their classical vertex correction). On the 
other hand, our detailed evaluation of the metric correction in gives us 
confidence in the correctness of our result. 



3.6 The vacuum polarization diagram 




(mi) iRa» 

ki i 


•Rsai (m 2 ) 



\ 


6 (a) 6(b) 

Figure 6: The vacuum polarization diagrams which contribute to the po¬ 
tential. Note that there exists a ghost diagram along with the graviton 
loop. 

The result for the vacuum polarization contribution can be found from 

is correctly pointed out in that the coefficient of the two-graviton vertex quoted 
in is too small by a factor of two. However, the numerical result for the loop 

integrals given therein is correct because of the presence of this symmetry factor. 
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the amplitude 


jppaXi 

Mfi(a)+6{b){q) = Tpa{ki,k2,mi)—^p—Ilx^f,u{q)—^p—r^5{k3,k4,m2) (38) 

where the vacuum polarization tensor is found from the effective Lagrangian 
obtained by’t Hooft and Veltman |3|, |^, — 


and is given by 


2G 

= -log(-^^) 

TT 


21 23 

''lafiV'rS 


23 2/ 

-^q {lapq-tqs + ri^5qaq0) 
21 


(40) 


11 


- Y4o^ iqaqsV0'r + q0qsva'y + qaq'rr]0s + q0q'rr]as) + ^qaq0q^qs 


(41) 


Contracting the various indices, we find a result 

43 

Me(a)+6{b) (^ = Y^G^"ii"i 2 log f (42) 

which is equivalent to that originally derived in |^, |^. After Fourier trans¬ 
forming we find a contribution to the scattering potential 


l^(a)+ 6 ( 6 ) (^) 


43 mim2G‘^ 
30 vrr^ 


(43) 


which is in agreement with that given by |15||. 


4 The result for the gravitational potential 


Adding up all the corrections to the non-relativistic potential we have our 
final result 


V{r) = - 


Gmim2 


G{mi +m 2 ) 41 Gh 

r IOtt 


( 44 ) 


The classical term in this potential agrees with Eq. 2.5 of IwasakifTT]]. 
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However, we note that the quantum component of the potential is not 
equivalent to any previous published result g H, 0, P, [I, H. As de¬ 
scribed above, we believe this result to be to definitive form of the non- 


relativistic scattering matrix potential. Our results agree with |^| for all 
diagrams except the vertex corrections^ However we have multiple checks 
on the vertex correction, as we have seen that our version leads exactly to 
the required form of the classical Schwarzschild metric, for both fermions 
and bosons. 


4.1 A Potential Ambiguity 

There is an important issue concerning the potential that has not been thus- 
far much discussed in the literature—that both the classical and quantum 
corrections to the potential are in some sense ambiguous. For the classical 
correction, this realization goes back to the papers of |^, 0, ^|, and these 
arguments generalize readily to the quantum component. In this section we 
discuss this ambiguity and argue that at one loop order our calculation of 
the quantum correction does have a well defined meaning. 

As explored in [0] the classical post-Newtonian potential is not invariant 
under a coordinate transformation of the form 


1 -|- a 


G{mi + m 2 ) 


(45) 


We are always free to make such a coordinate change, and that given in Eq. 
P modifies the classical correction to the potential via 


Gmim2 


l + c 


G{mi + m2) 


Gmim2 


1 -|- (c — a) 


G{mi + m 2 ) 


(46) 


Following 0, 0[ the ambiguity can be seen to arise in a field theory calcu¬ 
lation through a modification of the graviton propagator 


1 


1 


^0^ 




(47) 


accomplished by using the energy conservation in the following identical way 
which holds true for general x. 


1 


(1 - x)((p20 - Plo)'^ + (Pio - P3o)^)) - x{p20 - Plo){PiO “ P3o) “ f 


(48) 


®To be precise we are comparing to the third version of ||^, available in the electronic 
archive. We had several disagreements with the original version, all of which except the 
vertex correction have been corrected in the third version. 
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Using this propagator in order to derive the gravitational potential, the 
result will in general depend on x, with the relation to the coordinate change 
being a = — x). Not only will the terms of order ^ depend on x, 

but so will the corrections of order (c)' This indicates that the post- 
Newtonian part of the static gravitational potential is not well-defined in 
and of itself. 

The coordinate ambiguity also generalizes to the quantum part of the 
potential. There exists a coordinate redefinition 

(49) 


1 + f ^—2 


which changes the coefficient of the quantum term in the potential 


Gmim2 ,Gh 
- 1 -|- d—Y 


Gmim2 

r 


l + {d-P) 


Gh 

J.2 


(50) 


Therefore, we need to address the uniqueness of the quantum correction to 
the potential. 

Since general relativity is invariant under coordinate shifts, if the po¬ 
tentials are changed there must exist other modifications that compensate 
for these changes, leaving the resulting physics invariant. Within a Hamilto¬ 
nian formulation, these modifications take the form of momentum-dependent 
pieces—in order to make the potential well defined, one must also specify 
the momentum terms in the Hamiltonian. 

First we examine the classical ambiguity. Classically there are two di¬ 
mensionless variables available for the expansion away from the Newtonian 
limit 

Gm 
r 

In bound state problems these quantities appear at the same order, as can 
be seen by use of the virial theorem. In the Hamiltonian the leading terms 
are order 

Hi ~ —, - 

m r 

and are unambiguous. At next order we have pieces of the form 




p‘^ Gm? Gm 

n.2 ~ — o ,- 

m^ r r 


'^)(^) 

' r ^ m 


(53) 


There exists an ambiguity between the last two of these terms arising from a 
coordinate change, but this effect cancels for physical observables, as shown 
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explicitly in Writing the Hamiltonian in the center of mass frame to 
post-Newtonian order, we have 


H 


+ 


P 


{— ■ 

\2mi 2m2 

Gmim2 


1 “t" Ci- 


-(— + —) 

\ 8mf SmI J 

, (P-f)^ , G{mi +m2) 
+ 0 -h c- 


mim2 


mim2 


(54) 


In the standard Einstein-Infeld-Hoffmann coordinates, the coefficients a, b, c 
have the values 


a = 


b = 


1 + 3 


(mi + m 2 f 

mim2 


c = — 


(55) 


but under the coordinate transformation of Eq. ^ one finds the modifica¬ 
tions 


1 + (3 + 2a) 


(mi + m2)^ 


mim 2 


(mi + m 2 ) 

- a - 

2 mim2 

1 

-- -a 


(56) 


Therefore, in order to specify the correction to the static potential, one 
needs also to identify the momentum coordinates. In particular, the classi¬ 
cal potential calculated by Iwasaki is that appropriate for Einstein-Infeld- 
Hoffmann coordinates. 

This classical ambiguity has not been worked out explicitly at the follow¬ 
ing order in the expansion, but the general pattern is apparent. Specifically, 
at the next order one has four types of terms 


Hz ~ 


V 


5 



f Gm \ ^ Gm? / Gm \ ^ 

m \ r ) ' r \ r J 


(57) 


The last term here is noteworthy because it goes like 1/r^— i.e. like the 
quantum correction in the potential—but it is distinguishable by its mass 
dependence plus the fact that it is order G^. Such terms will also be ambigu¬ 
ous under a coordinate transformation, but the effects will cancel among the 
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effects terms of this order. So the “classical” coordinate change above does 
change the 1/r^ term in the potential, but has a specific form and the effects 
cancels among other classical terms in the Hamiltonian. 

Now let us examine the quantum effects, keeping only one power of the 
quantum expansion parameter 


Gh 

^2 

The quantum potential will enter the Hamiltonian at order 


(58) 


H,= 


Gm"^ Gh 


and the other term of this order is 


p‘^Gh 
~ m r2 


(59) 


(60) 


If one makes the “quantum” coordinate change of Eq. ^ , this will generate 
such a term in the Hamiltonian, which for physical observables cancels the 
effect of the change in the potential. Therefore in order to make the quan¬ 
tum correction to the potential well defined, one must specify the value of 
the terms contained in Hqp. However, the important point is that, in any 
coordinates in which one calculates quantum corrections, one will not gener¬ 
ate a term such as Hqp. This is because all quantum effects that involve the 
interparticle separation r arise from loop diagrams of order G^ while Hqp 
has only a single power of G. Quantum effects on a single particle line are 
of order G but do not involve r, while those involving two lines are of order 
G^. Since all calculations yield no term of the form Hqp, they must give the 
same quantum potential. The fully specified quantum potential is that one 
in the coordinates where Hqp vanishes. This, of course, is then the one that 
we have calculated. We conclude that, despite the expected general coordi¬ 
nate invariance, the calculated quantum potential is a well defined quantity 
to this order. 


5 Discussion 

The scattering amplitude has been used to provide a solid definition of the 
quantum corrections to the Newtonian potential. Our basic calculation is 
of certain nonanalytic terms in momentum space, the \/—q^ and (f' In —cf' 
terms. A transformation of these terms to coordinate space allows us to 
interpret these as long distance corrections to the potential. Our result is 
displayed in Eq. 
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The quantum corrections that we find for the scattering potential hold 
equally for the bound state potential. This is because the transformation 
between the two has no quantum component. The difference between the 
two is a classical correction that comes from iterating the lowest order po¬ 
tential. Dimensional analysis reveals that in the nonrelativistic limit this 
iteration has the ability to generate only a classical effect and not a quan¬ 
tum correction. 

We have found a result for the nonrelativistic potential which we believe 
is the final and complete result for this quantity. The potential matches 
the expectations from dimensional analysis as discussed previously i0] 
and the known ambiguity of the form of the classical correction has been 
seen to originate from the possibility of rewriting a potential energy term in 
the Hamiltonian in terms of kinetic energy and vice versa, as also discussed 
in 115, 19]. Such rewritings are not possible to do for the quantum term 

at the order to which we work. Therefore the quantum correction to the 
potential is a definite exact quantity. 

The quantum corrections are too small to be observed experimentally. 
However, the fact that these are reliably predicted is important for our 
understanding of quantum gravity. These effects are due to the low energy 
propagation of massless degrees of freedom and hence are uniquely predicted 
for any quantum theory of gravity that reduces to general relativity in the 
low energy limit^ In this sense, these are low energy theorems of quantum 
gravity. 


Acknowledgments 

N.E.J. Bjerrum-Bohr would like to thank the Department of Physics and 
Astronomy at UCLA for its kind hospitality and P.H. Damgaard for dis¬ 
cussions. The work of BRH an JFD is supported in part by the National 
Science Foundation under award PHY-98-01875. 

Appendices 


A Vertices and Propagators 

We begin by listing the Feynman rules which are employed in our calculation. 
For a derivation of these forms, see [^|. 

^Indeed, related effects have been found within the context of M-theory although 
the precise coefficients have not been calculated 
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A.l Scalar propagator 

The massive scalar propagator is: 



i 

q2 _ jyiP, _j_ 


A.2 Graviton propagator 

The graviton propagator in harmonic gauge can be written in the form; 


where 


a/3 ^5 

Q 


+ ie 


■paP^S _ 

2 




A.3 2-scalar-1-graviton vertex 

The 2-scalar-1-graviton vertex is discussed in the literature. We write it as: 


where 


p 

fiu / 

o 'n 


p 


= Trip,p',m) 


Ti'^{p,p', m) = —^ [p^p'’^ + p’^p'^ — 7]^" [{p ■ p') — m^)] 


A.4 2-scalar-2-graviton vertex 

The 2-scalar-2-graviton vertex is also discussed in the literature. We write 
it here with the full symmetry of the two gravitons; 



_ rjXpa 
— '2 


ip,p',m) 
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where 




_ _ ^^.XjpaaP ^ I + p'^p^) 


with 


_1 |^,Ap._l^,A^p.j^[(p.y)_^2] 


^a/S'fS niVct-yV/BS + ^aSV/S'f)- 


(61) 


A.5 3-graviton vertex 

The 3-graviton vertex can be derived via the background field method and 
has the formP, [ 1 l^] 





V 


where 

fiiy 


r3Z,sik,q)=-'-^x{r^ 


2 '' \ ' 0/375 


+ ‘^qxQ< 

+ 


+ {k- qY{k - qY + qf^q^ - ^r]^'^q‘^ 

T aX j I T crX T j ^cr j i/X j fia j uX 

^af) ^■yS ' ^'fS ~ ^a/3 ^'fS ~ ^yS ^afj 

qxq^ (^VapIyS + riyslap"^^ + q^q’' [r]al3Y 


'■yS ^ + V'ySlaP 

q^ (valBlyS - V'ySlap^"') “ 


+ [2gA(4/3""+ 4/"- 45 - lyS ""4/5 /^") 

+(4/3/45 + lo ^ rw ) +^" 4 .^a ( 4 /"a 
+ I (^^ + (^ “ q)‘^) [4/3^'^4 


r '^4-7' ^Pj '^1 

75p ' ^yS ^ci0p J 

^a/"45/ + 45 ^"4/5/ - 

^yS ^'"Vapk^ + lap^’^V'Tsik - qf) I) 


(62) 
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B Useful Integrals 

In the evaluation of the various diagrams we employ the following integrals 


J = 




Jfiu - 


as well as 
I = i 


dH 1 _ ^ r or] 

(27r)4 + J + 


dH 


1. 


(27r)4 P{1 + g)2 
dH 

(27r)4 P{1 + q)^ 

' dH 


327r2 L 
i 

327r2 




+ 


- -L) - q ri^A - -L 


(63) 

(64) 
+(.65) 


1 


(27r)^ P{1 + qY{{l + /c)2 — rm?) 

i 


I„ = 


327r2m2 

f dH 


L 


(27r)^ P{1 + qy{{l + /c)2 — m2) 


IflU - 


327r2m2 

f dH 


kn 


- 1 - 


1 

2 m 2 


L- 


1 

4 m2 


H + q^i' L + -S 


+ 




(27r)^ P{1 + qy{{l + /c)2 — m2) 


327r2m2 

+ 


q^qu [ - L - -S + - 


lil 

2 m 2 


■L- 


lil 

8 m2 


■5 


1 1 


3 q 


o + )-^ + ) +9 Q-^ 


16 m 


( 66 ) 


(67) 


+ 


( 68 ) 
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(27r)^ + qY{{l + — w?) 


327r^m^ 


q^Qi'Qa { 1 ^ k^ki/kct I 


“1“ {q^ki/kct -\- Qi/k^ko^ + Qak^kjy^ 


1 < 7 ^ 1 

3 m^ lb m^ 


+ {Qf^guka + Qt^gaku + gugakf,) 

(^^fii/ka ~\~ ^fiaku + rji/aku^ ( 12 ^ 

+ {vt^uga + Vfiagu + Vuagti) ( - ^g^L - ■^g^s'j + ■■■ 


where we have defined L = ln(—( 7 ^) and S = ^ ^ . Only the lowest order 


nonanalytic terms are included in the above forms. Higher order nonanalytic 

contributions as well as the neglected analytic terms are denoted by the 

2 

ellipses. The following identities can be verified on shell, k ■ g = where 

k — k' = g and k‘^ = m? = k'"^. In some cases the integrals are used with k 

2 

replaced by —k', where k' ■ g = 

In order to do evaluate the box diagrams the following lowest order 
integrals are needed. The higher order contributions of non-analytic terms 
as well as neglected analytic terms are again denoted by ellipses. 


K = 


K' = 


J (27r)4 P{1 + gy{{l + /ci)2 - mf)((/ - ks^ - mf) 

i A ^ _l_ 

167r^mim2(7^ V 3mim2/ 

r ciH _ 1 _ 

J (27r)^ P {1 + gY{{l + kiY - m\){{l + k^Y - ml) 

i . IT \ 1 


167r^mim2(7^ 


- 1 + 


3mim2 


L +... 


Here ki ■ g = k 2 ■ g = k^ ■ g = —^ and k^ - g = ^, where ki — k 2 = 
ki — ks = g and A:f = mf = k^ together with k^ = ml = k^. Also, we have 
defined w = {ki ■ k^) — mim 2 and W = {ki ■ k^) — mim 2 - 

The following constraints for the nonanalytic terms of the above integrals 
hold true on shell: 

= 0 (72) 
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,2 „2 


( 73 ) 



(74) 


References 

[1] S. Weinberg, “Phenomenological Lagrangians,” PhysicaA 96 (1979) 


327. 


[2] See, e.g., J. Gasser and H. Leutwyler, Ann. Phys. (NY) 158, 142 (1984); 
Nucl. Phys. B250, 465 (1985). 

[3] G. ’t Hooft and M. Veltman, 69-94, Ann. Inst. Henri Poincare, Sec. 
A20, 1974. 

[4] M. Veltman, Gravitation, 266-327, edited by R. Balian and J. Zinn- 
Justin, eds., Les Houches, Session XXVIII, 1975, North Holland Pub¬ 
lishing Company, 1976. 

[5] M. H. Goroff and A. Sagnotti, “The Ultraviolet Behavior Of Einstein 
Gravity,” Nucl. Phys. B 266 (1986) 709. 

[6] S. Deser and P. van Nieuwenhuizen, “One Loop Divergences Of Quan¬ 
tized Einstein-Maxwell Eields,” Phys. Rev. D 10 (1974) 401. 

[7] S. Deser and P. van Nieuwenhuizen, “Nonrenormalizability Of The 
Quantized Dirac - Einstein System,” Phys. Rev. D 10 (1974) 411. 

[8] B. S. Dewitt, “Quantum Theory Of Gravity. 1. The Canonical Theory,” 
“Quantum Theory Of Gravity. li. The Manifestly Covariant Theory,” 
“Quantum Theory Of Gravity. lii. Applications Of The Covariant The¬ 
ory,” Phys. Rev. 160 (1967), 1113; 162 (1967) 1195; 162 (1967) 1239. 

[9] J. E. Donoghue, “Leading quantum correction to the Newtonian poten¬ 
tial,” Phys. Rev. Lett. 72 (1994) 2996 [|arXiv:gr-qc/93106^] . 

[10] J. E. Donoghue, “General Relativity As An Effective Field Theory: The 
Leading Quantum Corrections,” Phys. Rev. D 50 (1994) 3874 ||arXiv:gr- 
qc/9405^ . 


24 







[11] I. J. Muzinich and S. Vokos, “Long range forces in quantum gravity,” 
Phys. Rev. D 52, 3472 (1995) |arXiv:hep-th/9501083i . 

[12] H. W. Hamber and S. Liu, “On the quantum corrections to the Newto¬ 
nian potential,” Phys. Lett. B 357, 51 (1995) [|arXiv:hep-th/9505182 l. 

[13] A. A. Akhundov, S. Bellucci and A. Shiekh, “Gravitational interaction 
to one loop in effective quantum gravity,” Phys. Lett. B 395, 16 (1997) 
| arXiv:gr-qc/9611018 l. 

[14] N. E. J. Bjerrum-Bohr, Cand. Scient. Thesis, Univ. of Copenhagen 

( 2001 ). 

[15] LB: Khriplovich and G.G. Kirilin, “Quantum power correction to the 
Newton law,” [|arXiv:gr-qc/0207fT^ 

[16] S.N. Gupta,“Quantization of Einstein’s Gravitational Field: General 
Treatment,” Proc. Phys. Soc. London A65, 608, 1952; S.N. Gupta and 
S.F. Radford, “Quantum field-theoretical electromagnetic and gravita¬ 
tional two- particle potentials,” Phys. Rev. D, Vol 21, No. 8, 2213, 
1980. 


[17] Y. Iwasaki, “Quantum Theory of Gravitation vs. Classical Theory,” 
Prog. Theor. Phys. 46, 1587 (1971). 

[18] K. Hilda and H. Okamura, “Gauge Transformation and Gravitational 
Potentials,” Prog. Theor. Phys. 47, 1743 (1972). 

[19] B. M. Barker and R. F. O’Connell, “Postnewtonian Two Body And N 
Body Problems With Electric Charge In General Relativity,” J. Math. 
Phys. 18, 1818 (1977) [Erratum-ibid. 19, 1231 (1978)]. 

[20] B. M. Barker and R. F. O’Connell, “Gravitational Two Body Problem 
With Arbitrary Masses, Spins, And Quadrupole Moments,” Phys. Rev. 
D 12 (1975) 329. 


[21] G. Modanese, “Potential energy in quantum gravity,” Nucl. Phys. B 
434 (1995) 697 []arXiv:hep-th/9408103] . 

[22] Diego A.R. Dalvit and Francisco D. Mazzitelli, Phys.Rev.D50:1001- 
1009,1994 |arXiv:gr-qc/9402003|] . 

[23] K. A. Kazakov, “On the notion of potential in quantum gravity,” Phys. 
Rev. D 63, 044004 (2001) |arXiv:hep-th/0009220i . 


25 


















[24] J. F. Donoghue, B. R. Holstein, B. Garbrecht and T. Konstandin, 
“Quantum corrections to the Reissner-Nordstroem and Kerr-Newman 
metrics,” Phys. Lett. B 529 (2002) 132 arXiv:hep-th/0112237 ]. 


[25] N. E. J. Bjerrum-Bohr, “Leading quantum gravitational corrections to 
scalar QED,” Phys. Rev. D 66, 084023 (2002) [ arXiv:hep-th/0206236 |. 


[26] N. E. J. Bjerrum-Bohr, J. E. Donoghue, and B. R. Holstein, ’’Quantum 
Corrections to the Schwartzschild and Kerr Metrics,” in preparation. 


[27] J. E. Donoghue and T. Torma, “On the power counting of loop dia¬ 
grams in general relativity,” Phys. Rev. D 54, 4963 (1996) [ arXiv:hep- 
I th/9602~T^. 


[28] K. Becker and M. Becker, “Quantum gravity corrections for 
Schwarzschild black holes,” Phys. Rev. D 60, 026003 (1999) [|arXiv:hep- 
th/98100^ . 


26 












